Model for curvature-driven pearling instability in membranes 
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A phase-field model for dealing with dynamic instabilities in membranes is presented. We use it 
to study curvature-driven pearling instability in vesicles induced by the anchorage of amphiphilic 
polymers on the membrane. Within this model, we obtain the morphological changes reported in 
recent experiments. The formation of a homogeneous pearled structure is achieved by consequent 
pearling of an initial cylindrical tube from the tip. For high enough concentration of anchors, 
we show theoretically that the homogeneous pearled shape is energetically less favorable than an 
inhomogeneous one, with a large sphere connected to an array of smaller spheres. 

PACS numbers: 87.16.Ac,87.16.Dg,87.68.+z,47.20.Ma 
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Introduction.- Intracellular transport in eukaryotic 
cells is an essential process in cell biology [if. Defor- 
mations of different intracellular membranes, as those 
from the Golgi apparatus or the endoplasmic reticulum, 
lead to budding of transport vesicles, which may eventu- 
ally undergo fission by the action of membrane proteins 
[3,Hl3- A membrane shape instability is, in some cases, 
the onset of those processes. Our work is thus motivated 
by these facts, aiming to give a simple physical interpre- 
tation of these morphodynamic effects in the cell interior. 

Dynamic instabilities in lipid vesicles and membranes 
have been widely studied experimentally. Tension-driven 
pearling induced by laser tweezers [11 , bilayer asymmetry 
p, polymer anchorage or osmotic perturbations [8| 
has been reported, as well as budding and tubulation 
induced by polymer anchorage j9|. 

A simple experimental model for membranes consists 
of a single component lipid bilayer and an amphiphilic 
polymer mimickingthe shaping effect of anchored mem- 
brane proteins 0, l^- When a certain amount of am- 
phiphilic polymer is introduced in the vicinity outside 
a lipid vesicle, its morphology changes destabilizing the 
equilibrium vesicle shape These experiments show 

that there is a coupling between polymer concentration 
on the membrane and local curvature. The polymer hy- 
drophobic backbones anchor to the outer leaflet of the bi- 
layer in order to minimize their hydrophobic interaction, 
acting thus as a wedge changing locally the curvature of 
the bilayer (see Fig. [1]). 

Theoretical works on tubular vesicle morphologies have 
been reported in the literature. Some of them have re- 
stricted themselves to find stationary shapes [lOl, [ll| . A 
hydrodynamical explanation for the pearling instability 
occurring when a sudden tension is induced in the mem- 
brane by laser tweezers was given in Refs. [H, [l^l- Some 
equilibrium models have dealt with shape instabilities 
due to the presence of anchoring molecules on the mem- 
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FIG. 1: Polymer wedge effect inducing a spontaneous cur- 
vature in a bilayer. A bilayer formed by one kind of lipids 
with zero spontaneous curvature tends to be flat (a). When 
a certain amount of anchor groups of an amphiphilic polymer 
gets stuck in the outer leaflet of the bilayer, a spontaneous 
curvature is induced (b). 



brane 0, [T^ . but a dynamic model has not yet been 
proposed. 

In this letter, we present a general method for deal- 
ing with morphological instabilities in membranes. In 
particular, we study the pearling instability induced by 
the anchorage of amphiphilic polymers on a lipid vesicle, 
in order to find and characterize the shape instabilities 
observed experimentally Q. We also show, motivated 
by our simulations, that the inhomogeneity in the pearl 
size which appears at high concentration of anchors on 
the membrane can be explained by energy considerations, 
consistently with the picture presented in [3]. 

Theoretical treatment. - A lipid vesicle subject to an ef- 
fect which causes an asymmetry between the two leaflets 
of the bilayer vesicle can be elastically described by the 
Canham-Helfrich Hamiltonian [l^, [l^ , 



Uc-n = |^(2H-Co)'d5, 



(1) 
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where k is the bending modulus, H is the mean curvature, 
and F is the vesicle surface. The asymmetry is introduced 
through the spontaneous curvature, Cq. Minimization of 
this Hamiltonian under the constraints of constant area 
and volume yields the equilibrium shapes of vesicles p7| . 
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To incorporate the effect of curvature generation by the 
anchored polymers, we assume a hnear coupHng between 
the spontaneous curvature and the polymer concentra- 
-(0) , r(i) ^(^ ^^ r(o) 



tion [IJ, Co(x,t) = C^"^ + C^^^p(x,t), where Q"^ is the 



bare spontaneous curvature, i.e. due to the asymmetry 
between the two leaflets of the bilayer before the poly- 
mer anchorage, Cq^'' is the polymer-induced spontaneous 
curvature, and p(x, t) is the local density of polymer. It 
has to be noted that here, the spontaneous curvature is, 
a priori, a dynamic non-homogeneous function. 

Phase-field model. - The dynamics of a membrane con- 
sidered as a surface is a paradigm of a moving boundary 
problem. A useful method for dealing with such prob- 
lems is the phase-field approach (see [19| for a review on 
surface tension-driven models). A phase-field model for 
taking into account the bending energy of lipid vesicles 
was presented in Refs. [l^, [2l|. The free energy func- 
tional proposed there was shown to be consistent with 
the already known stationary shapes of fluid vesicles. 

We present here an extension of the phase-field model 
for the Canham-Helfrich bending energy derived in 
Ref. [10], in order to take into account a polymer- 
dependent spontaneous curvature. The free energy 
^[4), p], is a functional of the local polymer density, p(x), 
and of a phase-field, (^(x, t), whose level-set {4>{x) = 0} 
locates the interface. 
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where 



e2v20(a;), (3) 



where Co is related to the spontaneous curvature as 
V^Co — Co, and may depend in principle on any dy- 
namic variable of the problem. Phase-field models are 
also called diffuse interface models, since the phase-field 
is a continuous albeit abruptly-changing function at the 
interface. The width of this diffuse interface is related to 
the small parameter e. In the sharp-interface limit, e — s- 
0, the phase-field model is equivalent to the equilibrium 
Canham-Helfrich Hamiltonian Eq. [1] As in Ref. [20| . 
we write an effective free energy functional JFcff [(/), p] = 
!F[cf), p] + Jj-j a[x)a[(f)]Ax, where a[(f)] — e | V(/)| , is the 

surface area density, and (t(x) a local Lagrange multiplier 
which ensures local area conservation. The conservation 
of the inner volume of the vesicle is achieved via a relax- 
ational conserved dynamics of the form ^ = i^^^f'^ > 
leading to the dynamic equation for the phase-field, 

^ = 2V2{(302-l + 2eCo(p)0)$[0]-e2v2$[</>] 



'a{x)V^4], 



(4) 



where (j(x) = ^a{x). A dynamic equation for this La- 
grange multiplier is also needed. A first-order Lagrangian 



method |22| is used to implement the conservation of the 
surface area. 

A dynamic equation for the density field, p(x, t), could 
also be straightforwardly introduced. However, for sake 
of simplicity, we are concerned here with the case of a 
homogeneous polymer distribution along the membrane, 
due to e.g. a global application of the polymer. A sys- 
tematic study for non-homogeneous polymer concentra- 
tion and polymer diffusion along the membrane si beyond 
the scope of this letter [SsJ . 

Eq. ([¥]) is a highly non-linear dynamic equation which 
must be solved numerically. We used a standard finite- 
difference scheme for the spatial discretization and an 
Eulcr method for the time-derivatives [20] . In all the re- 
sults shown in this letter we used the value of the small 
parameter, e, equal to the mesh size of the lattice. The re- 
sults are robust under variations of this parameter. Both 
the time and space discretizations are chosen in order 
to satisfy the Courant-Friedrichs-Lewy stability criterion 

Results and discussion.- In the experiments by Tsafrir 
et al. [7|], amphiphilic polymers are introduced in the 
bulk outside the vesicle, both globally and locally close 
to the tip of the tube. These molecules diffuse in the 
bulk until they come across the membrane, where they 
get stuck in such a way that their hydrophobic part an- 
chors in the bilayer (in order to satisfy the hydrophobic 
interaction). Once a polymer is anchored in the bilayer, 
it diffuses along the membrane. 

Here we consider the situation of global application of 
the polymer. In this case, we assume that the polymer 
concentration reaches a homogeneous stationary profile 
along the membrane almost immediately. The dynamic 
evolution is thus fully understood from the shape dy- 
namics, so there is no need for a dynamic equation for 
the density field. 

Onset of the instability.- Tsafrir et al. Q studied the 
pearling instability in tubes whose length is much larger 
than their diameter. Within our model, we can find the 
shape of the tube at the onset of the instability (see 
Fig. [21). The experimental and predicted shapes are in 
good qualitative and quantitative agreement with each 
other. For instance, we can measure the ratio between 
the radius of the first pearl and the radius of the neck 
connecting it with the tube, and see that it gives a value 
of about 3 in both the experiment and the simulation. 
We can thus assert that there is no need for an inho- 
mogeneous polymer concentration on the membrane in 
order to start the instability, but it can just be triggered 
by the global change of the preferred curvature. 

Low polymer concentration. - Further addition of poly- 
mer solution increases the concentration of anchor chains 
in the membrane, therefore increasing the induced spon- 
taneous curvature. Fig. [3] shows the time evolution of 
a semi-infinite cylindrical tube with an endcap and the 
other end connected to a lipid reservoir. The homoge- 
neous spontaneous curvature induced by the anchorage 
of the polymers made the initial cylindrical shape to be 



(a)Experimcntal Result (from Rof. 7]) 




(b)Phase-Field Simulation 

FIG. 2: Onset of the pearling instability. Comparison of the 
experimental result from Ref. JJ (a), and the phase-field simu- 
lation (b) . It is important to note that in the simulation there 
is no fitting parameter, but we just let the system evolve from 
an initial tubular shape, under a relatively low homogeneous 
induced spontaneous curvature, Co = 0.48, below the pearling 
instability limit. 
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(a) Experimental Result (from Ref. Q) 



(b)Phase-Field Simulation 



unstable and to create pearls. 

We define the volume to area ratio, A = V/A. For 
spontaneous curvatures between Co = 1/(2 A) and Co — 
2/(3 A), Deuling and Helfrich [U showed that there ex- 
ist minimal surfaces, called Delaunay surfaces, which are 
global minima of the bending energy Eq. ^ , for cylindri- 
cal shapes. The two limiting Delaunay shapes are a cylin- 
der (for Co = 1/(2 A)) and a set of spheres (Co — 2/ (3 A)). 
Unduloids, the one-parameter family of Delaunay shapes, 
interpolate smoothly between them. For higher spon- 
taneous curvatures, the condition of vanishing bending 
energy cannot be fulfilled, and therefore the stationary 
shape corresponds to a non-vanishing minimum of the 
free energy. Our simulations show that, when the poly- 
mer solution is applied globally, i.e. when there is a simul- 
taneous anchoring of polymer everywhere on the surface, 
the way of reaching a pearled structure from a cylindri- 
cal one is not by going through the family of equilibrium 
Delaunay shapes, but by a completely different dynamic 
evolution: creating pearls one by one from the tip of the 
tube (see Fig.[3i;b)). 

We checked how the free energy Eq. ^ varies in time 
seeing how each pearl formation is associated with an en- 
ergy barrier which is crossed by thermal activation (due 
to numerical noise in our simulations). It is important to 
remark that there occurs no fission in the tube, but there 
is a narrow neck joining any pair of pearls, as seen in the 
experiments. The width of this neck is of the order of the 
mesh size. There is therefore no change in the topology 
of the tube. The tip of the semi-infinite tube is the point 
where pearls start when the polymer solution is applied 
both locally on the tip and globally 

High polymer concentration. - It is seen in the exper- 
iments (see Fig. IH^a)) that when the concentration of 
anchored molecules on the membrane increases, there ap- 
pears a gradient in the size of the pearls. The higher the 
concentration, the more pronounced the size gradient is. 
This is so because the higher the amount of anchored 
polymer on the membrane, the higher the induced spon- 



FIG. 3: Dynamic evolution (time goes by from upper im- 
ages to down) of a tube with an polymer concentration on 
the membrane such that Co = 0.68. We show (a) the ex- 
perimental results from 0] and (b) the phase-field simulation 
for comparison. Pearls start at points where the tube loses 
its perfect cylindrical geometry, namely the cap. The simula- 
tion has been performed on a 200 x 40 axisymmetric lattice. 
No-flux boundary conditions at the lateral walls have been 
implemented. 



tancous curvature. As we mentioned before, there are no 
energy- vanishing surfaces for Co > 2/3. Tsafrir et al. [3| 
suggested that a shape consisting of a chain of equally 
sized Helfrich spheres connected to a larger sphere is, 
in terms of free energy, favorable to a chain of equally- 
sized spheres. Some pearls having the mean curvature 
equal to the spontaneous curvature are formed, but not 
all of them can fulfill this condition, under the volume 
and area constraints. Then, a larger sphere should form 
in order to keep these constraints. They argued that 
an inhomogeneity in the polymer concentration on the 
membrane (and thus in the distribution of the induced 
spontaneous curvature), may reduce the Helmholtz free 
energy F = 7ic-H — TS, T being the temperature and S 
the entropy. 

However, we found that for higher values of the spon- 
taneous curvature, a situation with a pearl size gradi- 
ent may be energetically favorable to that of a chain of 
equally-sized pearls, even with the same spontaneous cur- 
vature all along the vesicle (Fig. 2]) . It is interesting to 
note that for values of the spontaneous curvature just 
slightly higher than 2/3, the homogeneous case keeps its 
stability against the inhomogeneous case. Once a criti- 
cal value Cg'^ is achieved, the homogeneous case destabi- 
lizes against the inhomogeneous one (see Fig.|4l^c)). This 
critical value depends on the area of the vesicle or, in 
other words, on the length of the initial tube. This criti- 
cal value decreases with increasing tube length, reaching 
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FIG. 4: Inhomogeneous pearling, (a) Experimental result 
from Ref. 01, and (b) phase- field simulation using Co = 1.1. 
(c) Plot of the energy difference AE between the bending 
energy corresponding to a set of spheres, and the one associ- 
ated with a set of small spheres added to a bigger one, with 
respect to the increment of spontaneous curvature from the 
value Co — 2/3, when equally-sized spheres have zero bend- 
ing energy. For spontaneous curvatures bigger than a critical 
value, the homogeneous configuration is energetically less fa- 
vorable than the inhomogeneous one. The area of the vesicle 
is finite, and chosen in such a way that the homogeneous 
pearled chain consists of 8 spheres. Magnitudes are measured 
in normalized units, given by k = 1 and A = 1. 



2/3 in the case of infinite tubes. Therefore, we state that 
the inhomogeneous pearl size experimentally found by 
Tsafrir et al. , is not due to a inhomogeneous polymer 
distribution, but it is of a purely energetic nature. 



Conclusions." In summary, a phase-field model for 
dealing with dynamic instabilities in vesicles has been 
studied in the case of curvature-driven pearling instabil- 
ity induced by the anchorage of amphiphilic polymers on 
the bilayer. We showed that the morphological changes 
reported in the experiments are explained by the gen- 
eration of curvature by the anchors. We considered the 
situation in which polymer is applied globally, account- 
ing for the main experimental results. Calculations in the 
framework of the Canham-Helfrich model showed that for 
a high enough homogeneous concentration of anchors, 
the homogeneous pearled shape is energetically less fa- 
vorable than an inhomogeneous one, with a large sphere 
connected to an array of smaller spheres. 
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